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Title: Doppler Spread Estimation 
Field of the Invention 

The present invention relates to the estimation of Doppler spread and mobile 
5 stations speed. 

Background to the Invention 

A radio signal from a transmitter to a receiver will travel along different propagation 
paths as the signal is scattered by obstacles, such as houses and other objects. This 
10 leads to signals received with different time delays, so called multipath propagation. 
When the receiver starts moving Doppler shifts will be introduced. When several 
paths arrive at the same time delay a Doppler spread, or Doppler spectrum, is 
generated, 

IS Knowledge of the Doppler spread has been found to be useful for enhancing the 
operation of receivers as described, for example, in Morelli, M. a/,, "Further 
Results in Carrier Frequency Estimation for Transmissions Over Flat Fading 
Channels", IEEE Communications Letters vol. 2, no. 12, December 1998. The 
Doppler spread can also be used to infer the scalar speed of a receiver relative to 

20 local radio reflectors. This scalar speed is not the radial velocity that one can 

determine by measuring a Doppler shift. In most cases, the speed will correspond 
to the speed of the receiver, such as in the case of a receiver in a moving car 
surrounded by stationary street furniture. 

25 Summary of the Invention 

The present inventors have established that the Doppler spread associated with a 
transmission path can be estimated reliably from the derivative of the transmission 
path's transfer function. 

30 According to the present invention, there is provided a method of estimating the 
Doppler spread of a radio signal, the method comprising 
receiving a radio signal; 



T 



deriving a value for the derivative of the envelope of the path transfer 
function for said radio signal; and 

computing an estimate of the Doppler spread of said radio signal from said 
derivative value. 

5 

The transfer function of a path from a transmitter to a moving receiver is subject to 
three amplitude effects, namely path loss, slow fading and fast fading. It is the fast 
fading which is associated with Doppler spread. The fast fading operates at time 
scales many orders of magnitude smaller than path loss and slow fading. 

70 Accordingly, for a signal transmitted at a constant or slowly varying power, the 

derivative of the envelope of the received signal will be determined predominantly 
by the fast fading, i.e. Doppler spread, and the path loss and slow fading effects can 
be disregarded. Modulation effects can be compensated for by demodulating the 
received signal with a signal corresponding to the modulating signal which can be 

IS predetermined, for example by means of reference codes transmitted at 



predetermined times. 



However, the path transfer function can be estimated by techniques, other than by 
demodulating the received signal with the modulating signal, such as that described 
20 in Wu, J. a/., "Blind Channel Estimation Based on Subspace for Multicarrier 

CDMA", Vehicular Technology Conference, 2001. VTC 2001 Spring. IEEE VTS 
53rd , Volume: 4 , 2001. 

In the minimal case of an unmodulated carrier, the received signal's envelope will be 
2S proportional to the transfer function of the path. 

Preferably, therefore, said value for the derivative of said envelope is derived by 
low-pass filtering an envelope signal representing the path transfer function 
envelope to band limit it and filtering the band-limited envelope signal using an FIR 
30 filter. More preferably, said envelope signal comprises a sequence of samples 
representmg the path transfer function envelope. 
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Preferably, the computing of said estimate of the Doppler spread comprises 
determining the variance of said derivative value. More preferably, the computing 
of said estimate of the Doppler spread comprises determining a value indicative of 
the received power of said radio signal. Still more preferably, the Doppler spread 
estimate is calculated by determining the square root of the result of dividing twice 
said variance by said value indicative of the received power of the radio signal. 

Preferably, the Doppler spread estimate is calculated in accordance with the 
formula:* 



10 Doppler spread oc 

where 

and 



1 N 

15 where ris the magnitude of the radio signal. 
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According to the present invention, a Doppler spread, determined according to the 
present invention, is used to estimate the speed of a mobile station in a wireless 
communication system. The speed estimation method may perform processing 
20 applicable to different speed ranges, e.g a low speed range and a high speed range, 
in parallel and then select the more appropriate speed value as the final speed 
estimate. 

According to the present invention, there are also provided mobile stations 
25 including processing means configured, e.g. by software, for performing methods 
according to the present invention. 



Brief Description of the Drawings 

Figure 1 is a block diagram of a mobile station; 
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Figure 2 is a flowchart illustrating a speed estimating program according to the 
present invention; 

Figure 3 is a flowchart illustrating a signal conditioning process of the program of 
Figure 2; 

5 Figure 4 is a flowchart illustrating a speed estimating process of the program of 
Figure 2; and 

Figure 5 is a flowchart illustrating a speed estimate selection process of the program 
of Figure 2. 

Detailed Description of the Preferred Embodiment 

A preferred embodiment of the present invention will now be described, by way of 
example, with reference to the accompanying drawings. 

Referring to Figure 1, a WCDMA mobile station comprises an antenna 1, an rf 
subsystem 2, a baseband DSP (digital signal processor) subsystem 3, an analogue 
audio subsystem 4, a loudspeaker 5, a microphone 6, a controller 7, a liquid crystal 
display 8, a keypad 9, memory 10, a battery 11 and a power supply circuit 12. 

The rf subsystem 2 contains if and rf circuits of the mobile telephone's transmitter 
and receiver and a frequency synthesizer for tuning the mobile station's transmitter 
and receiver. The antenna 1 is coupled to the rf subsystem 2 for the reception and 
transmission of radio waves. 

The baseband DSP subsystem 3 is coupled to the rf subsystem 2 to receive 
25 baseband signals therefrom and for sending baseband modulation signals thereto. 

The baseband DSP subsystems 3 includes codec and RAKE functionaUty, which are 
well-known in the art, and is programmed for estimating a received signal Doppler 
spread and the speed of the mobile station. 

30 The analogue audio subsystem 4 is coupled to the baseband DSP subsystem 3 and 
receives demodulated audio therefrom. The analogue audio subsystem 4 ampUfies 
the demodulated audio and appUes it to the loudspeaker 5. Acoustic signals. 
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detected by the mictophone 6, are pre-amplified by the analogue audio subsystem 4 
and sent to the baseband DSP subsystem 3 for coding. 

The controller 7 controls the operation of the mobile telephone. It is coupled to 
S the rf subsystem 2 for supplying tuning instructions to the frequency synthesizer 
and to the baseband DSP subsystem 3 for supplying control data and management 
data for transmission. The controller 7 operates according to a program stored in 
the memory 10. The memory 10 is shown separately from the controller 7. 
However, it may be integrated with the controller 7. 

W 

The display device 8 is connected to the controller 7 for receiving control data and 
the keypad 9 is connected to the controller 7 for supplying user input data signals 
thereto. 

fS The battery 1 is connected to the power supply circuit 12 which provides regulated 
power at the various voltages used by the components of the mobile telephone. 

The controller 7 is programmed to control the mobile station for speech and data 
communication and with application programs, e.g. a WAP browser, which make 
20 use of the mobile station's data communication capabiHties. 

In order to estimate the speed of the mobile station, the DSP subsystem 3 is 
programmed to perform calculations appropriate for low and high speeds in parallel 
and then select the appropriate result. 

25 

Referring to Figure 2, the speed estimation program of the DSP subsystem 3 
receives a stream of in-phase and quadrature baseband signal components and 

provides samples time-domain of the transfer function (H) of the signal paths 
having the various delays handled by the RAKE processing. In the present 
30 example, these transfer function samples are obtained by demodulating the received 
signal for reference symbols with values representing the transmitted signal. 



Signal conditioning is petformed on the ttansfer function samples (step pi a and 
plb) using algorithms adapted for low and high speed ranges respectively. The low 
and high speed ranges overlap in the present example. However, the upper edge of 
the lower range may simply meet the lower edge of the upper range. The 
conditioned signals are then used to produce speed estimates (steps p2a and p2b) 
and the appropriate speed estimate is selected (step p3). 

The signal conditioning processes pi a, plb loop through the RAKE fingers which 

process transfer function signals for respective different path delays. Of course, if 
the mobile station does not employ a RAKE system, such as where the mobile 
station is not a CDMA device, there will be no need for this loop. 

Referring to Figure 3, in the low speed signal conditioning estimation process pi a, it 
is first determined whether the RAKE finger to be processed is locked (step si). A 
"locked" finger is one whose allocated Ufetime has expired. If the finger is locked 
at step si, processing for the current finger is skipped. 

However, if the current finger is not locked at step si, the state of the finger, i.e. 
deallocated or allocated to a radio path, is determined (step s2). 

In the case of the current finger being in the deallocated state at step s2, it is 
determined whether the finger's lifetime has expired (step s3). If the finger's 
hfetime has expired at step s3, the finger is locked (step s4) and further processing 
for It is skipped. However, if the finger's lifetime has not expired at step s3, all 
variables and parameters for the current finger are reset (step s5). 

In the case of the current finger being allocated at step s2 and following step s5, it is 
determined whether a valid reference symbol has been detected (step s6). The 
reference symbols occur regularly in signals transmitted to the mobile station and 
have known and values. Therefore, values associated with like signals are used for 
calculation of the Doppler spread and speed. In the present case, the power at 
which the reference signals are transmitted is constant. However, if the power at 
which they are transmitted varies sufficiently quickly to affect the determined 



- 7 - 

derivative value for the envelope, the magnitude of the envelope can be normalised 
using information about their transmission power. This information may be 
provided in control channels, for example. 

5 If a valid reference symbol has been detected, in-phase (I) and quadrature (Q) 
component values for the transfer function are obtained (step s7). When the 
aforementioned I and Q values have been obtained, they are used to calculated the 
magnitude of the transfer function's envelope (r) (step s8). 

10 However, if a valid reference symbol has not been detected, the previous value of 
transfer function's envelope magnitude is used to avoid gaps in the sequence of 
Zl envelope magnitude values (step s9). 

g/j When the envelope magnitude value has been obtained, it is determined whether 

^9 15 continuous invocation mode is being used (step slO). Continuous invocation mode 
is the normal mode of operation with an envelope value be obtained for each slot, 
i.e. regularly and frequently. However, under some circumstances, the speed 
fU estimation program may not be run in some slots. For instance the program may 

V- not be called during the transmission gap in compressed mode. 

H 

If operation is continuous, the new envelope magnitude value is fed into a low-pass 
filtering process (step sll). The low-pass filtering is provided by implementing a 3'"* 
order Butterworth IIR filter with a cutoff at SOOHz. The filtering process for the 
low speed thread also performs downsampling by a factor of 2. The bandwidth of 
25 the filter is set to twice the maximum Doppler shift that can occur in the speed 

range covered. Thus, SOOHz corresponds approximately to a speed of 125km/h at 
2.17GHz. 

If the operation is discontinuous, the filtering process sll is skipped and the output 
30 of the filtering process sll is replaced with the current envelope magnitude value. 



The result of the filtering process sll or the current envelope magnitude, as the 
case may be, is used to calculate the square of the envelope magnitude (r/), the 



X 



derivative of the envelope magnitude (rj and the square of the derivative of the 
envelope magnitude {r^)- Since, the envelope magnitudes are time-spaced samples, 
the derivative is approximated in the low speed thread using: - 

f- = ^< ~ ^-Ar 

^T 

5 where ^ indicates an approximate or estimated value and AT is the sampUng 
interval. It can be seen that this is a two-tap FIR filter. 

These values are used to update accumulated values ^^^nT'>^^KT->^^^nT ^he 

n=l 71=1 n=\ 

square of the envelope magnitude, the approximate derivative of the envelope 
10 magnitude and the square of the derivative of the envelope magnitude, where N is 
n an observation window length in terms of sample. In the present example, N is 

R 750. 

s 

f5 If all of the fingers have been processed at step sl2, the signal conditioning is 

^ 15 complete. Otherwise, the process returns to step si. 

Vz The high speed signal conditioning plb is the same as the low speed signal 

4- conditioning except that the low-pass filter has a cut-off at IkHz rather than 

P 

fy SOOkHz and the envelope magnitude derivative r, is approximated using:- 

which is a 3 -tap FIR filter having a zero coefficient for the middle tap. 

The bandwidth of the filter is set to twice the maximum Doppler shift that can 
occur in the speed range covered. Thus, IkHz corresponds approximately to a 
25 speed of 250km/h at 2.17GHz. 

The results of the low speed and high speed signal conditioning process pla, plb 
are supplied to the low speed and high speed speed estimating processes p2a, p2b 
respectively. 

30 



Referring to Figure 4, at the start of the low speed speed estimating process p2a, an 
initialisation routine (step s21) sets a "finger valid flag" to true and a "maximum 
energy" variable to —1. After initiahsation, the process enters a loop so that all of 
the RAKE fingers are processed in turn. 

At the start of the loop, it is determined whether the current finger's Ufetime is 
greater than or equal to a threshold (step s2). The threshold is set to identify 
fingers what are sufficiently old to provide sufficiently good data for reliable speed 
estimation. If the current finger's Ufetime is less than the threshold, the process 
moves on to the next finger. 

If at step s22, the current finger is sufficiently mature, a value Z?q and an energy 

related value for the current finger are calculated (step s23). is indicative of the 
power of the current finger and is calculated using:- 

N 

where E{r^} is the mean of ^^r^f calculated at step sll in Figure 3. 

The energy related value for the current finger is calculated by multiplying Bq by the 
current finger's age. 

Following step s23, it is determined whether the current finger has the largest 
energy related value of the fingers processed so far (step s24) by comparing its 
energy with the "maximum energy" variable. If this is not the case, the process 
moves on to the next finger. However, if this is the case, the "finger vahd flag" is 
set to true to indicate that data is now available for speed estimation, the "maximum 
energy" variable is set to the energy of the current finger and an index, identifying 
the current finger, is stored (step s25). 

When all of the fingers have been processed (step s26), it is determined whether 
data suitable for speed estimation has been obtained by checking the state of the 
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"finger valid flag" (step s27). If the "finger valid flag" is false, no further 
calculations are performed and the result of the low speed speed estimating process 
p2a is noted as being unrehable by setting a reliability flag to false. However, if the 

"finger valid flag" is true, the Doppler spread f^^p.^^ estimated for the last finger 

whose index was stored at step s25 using:- 




where 

b,=E{}^}-{E{i'}f 

N ^ 

where E{r^} is the mean of ^r^j- calculated at step sll of Figure 3 and E{r} is the 
mean of ^r^y- calculated at step sll of Figure 3 (step s29). 

The newly calculated Doppler spread estimate is checked to determine whether it 
falls within a range of valid values (step s30). If the Doppler spread estimate is 
outside the vahd values range, the process proceeds to step s28. However, if it is 
within the valid values range, it is used to computer a speed estimate for the mobile 
station using: - 




where c is the velocity of light in free space and is the carrier frequency of the 
signal being received. It is determined whether the newly estimated speed value is 
within a valid range, i.e. below 500km/h in the present example, (step s32). If not, 
the process moves to step s28. If the speed value is valid, the reUability flag is set 
to true. 

The high speed speed estimating process p2b is identical except for the vaHd value 
ranges used at steps s30 and s32. 



Referring to Figure 5, when the speed estimation processes p2a, p2b have been 
completed, the speed selection process p3 determines whether both speed estimate 
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process have produced reliable estimates (step s41). If one or both have not 
produced reliable estimates, as indicated by their respective reUability flags, no 
speed estimate is selected, a speed estimate reliabihty flag is set to false (step s42) 
and the process terminates. 

If both speed estimation processes p2a, p2b produced valid speed estimates, the 
estimates are compensated for systematic errors introduced by the low-pass filtering 
and the speed estimation algorithm for low speeds (step s43). The compensation 
algorithm is a standard linear mapping of the form:- 

and in the present example a=1.7 and p=-12 for the low speed speed estimate and 
a=l.l and P— -20 for the high speed speed estimate. 

It is then determined whether the compensated low speed speed estimate is less 
than a threshold located within the overlap between the low and high speed ranges 
(step s44). If the low speed speed estimate is less than the threshold, it is 
determined whether the estimate is less than 0 (step s45) and, if so, the final speed 
estimate is set to 0 (step s46). If the compensated low speed speed estimate is not 
less than 0 at step s45, the final speed estimate is set to the compensated low speed 
speed estimate (step s47). 

If the compensated low speed speed estimate is not below the threshold at step s44, 
the final speed estimate is set to the compensated high speed speed estimate (step 
s48). 

Following steps s46, s47 and s48, the speed estimate reUabiUty flag is set to true. 

If the speed estimate reUabihty flag is true, a speed estimate will now be available 
for use by the mobile station and can be reported to a fixed network node where it 
may be used to control handovers, for example. 

It will be appreciated that many modifications may be made to the preferred 
embodiment described above. For example, additional loops may be introduced 
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into the signal conditioning and speed estimating processes for iterating through 
different signal sources in the case of transmit diversity being used. 

In the foregoing, speed estimates are determined for two speed ranges. However, 
5 calculations may be performed for more or fewer speed ranges as required. More 
ranges give better performance but at the cost of increased complexity. Also, 
higher order low-pass filters give better performance but again at the cost of 
increased complexity. 

10 In the foregoing, the Doppler spread and speed calculations are performed in a 

mobile station. However, these calculations could be performed at a fixed network 
node using signal strength reports from a mobile station. Alternatively, the 
calculations could be performed in part by a mobile station, e.g. the signal 
conditioning processes, and in part by a fixed network node, e.g. the speed 

15 estimation and selection processes. 



